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ABSTRACT 



A single Fourier transform optical system was constructed and tested 
on two projects. First, film grain in a single Fourier transform operation 
was studied and a mathematical model of the film grain and a semi-quan- 
titative prediction of the grain's transform were developed and verified 
using the experimental system. The film grain was found to have a noise- 
like spatial frequency spectrum which was a decreasing function of fre- 
quency, and that this spectrum could be considered Gaussian only in 
the case of fine grain films containing binary signal images. Second, 
the system was used to investigate. a proposed method of detecting a 
ship wave. A deterministic relation between a ship wave and the Fourier 
transform of its photograph was shown to exist, A sample photograph of 
a ship wave was inserted into the system and operated on by a matched 
filter rotated in a plane normal to the optical axis , The output light 
intensity was determined as a function of the filter angular position. 

The output intensity when considered to be periodic was found to have a 
discernible spectral content due to the ship wave. 
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I. INTRODUCTION 



A. BACKGROUND 

Optical systems have received increased attention from physicists 
and engineers. The reason lies in the Fourier transform properties of a 
coherent optical system, and the laser has recently provided a practical 
source of coherent light. The optical image possesses two degrees of free- 
dom (i.e.; the coordinates of a point on the surface) in contrast to the 
single degree of freedom, time, as found in electrical systems. With 
two degrees of freedom, an optical system appears to offer distinct advan- 
tages over an equivalent electrical system in some applications . 

Since the electrical engineer is generally not familar with the 
Fourier properties of an optical lens, a brief description is presented 
here. More detailed analyses are given in Preston [l ], Francon [2], 
and Goodman [3]. Consider the system diagrammed in Figure 1. The 
input plane, PI, is illuminated by a planar, unpolarized, unmodulated, 
monochromatic light wave described by the equation: 

fjd) (1-1) 

where tc is the radian frequency of the light wave. The plane, PI, 
contains a light modulator t(x^,y 3 ^) such as a photographic transparency. 

The transmitted light is: 

/t) = fot (xi,yi)e^^’^ (1-2) 

where t is complex. Equation (1-2) may be written as: 

fi(xi,yi,t) = fot(x^,yJexp[ja.'t+ (|) (X3^,y^)] (1-3) 
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Assume the lens to be ideal and the aperture, la.'ge so that vignetting 

does not occur. Geometric optics shows that the energy diffracted at PI 

in the direction, v, will be summed at a point in the back, focal plane, P2 . 

It may be shown that the phase at the summing point, {y^z,Yz) , is: 

(I* (xg/ys) = 4(xi,yi)+ ^ (xiXg+yiys) + (1-4) 

where Xis the wavelength of light, R is shown in Figure 1 , and d is the 

optical path length from the origin of PI to the point (x 2 ,y 2 ). The light 

incident at the point, {^ 2 >y^ > becomes: 

12(^2 /y 2 /t)=foS^^'Vj' t(xi,yi)exp j[(J(xi,yi) + 

PI 

(x^xs+yiys) + ] dxjdyi (1-5) 

The phase modulation term, (| is usually made constant in prac- 

tice and most often detection is by energy so the term outside the integral 
may be replaced by its rms value. For small diffraction angles where R 

is approximately equal to the focal length, F, equation (1-5) reduces to: 

2 7T 

f2(x2/y2)=fa JJ t(x^,yjexp j — (xiX2+yiy2)dXi,dyi ( 1 - 6 ) 

PI 

where f„=^^exp j [^(| (xi,yi)+ 



Equation (1-6) is the two-dimensional, bounded, Fourier transform of the 
function, t(xj^y^) , multiplied by a system constant. This may be more 
readily apparent by rewriting the kernel function exp j (27 t/X F) (x^^ X 2 +yiy 2 ) 
as exp -j (cOj^x^+cOyyx ) and defining: 



OJ 



X 



60 = - 
y 



2tt 

77T Xp 

XF ^ 
2tt 



XF 



y2 



(1-7) 
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Thus a' snd are spacial frequencies defined in the back focal plane, 

X y 

P2 . Two successive optical transform operations will map the input image 
into the output with only a coordinate reversal and a constant amplitude 
change . 

B. APPLICATIONS 

One of the more obvious applications follows directly from the above. 
Since the plane of the first transform, or frequency, plane displays the 
spatial frequency content of the input function, t(xj^,yi), frequency 
modifying operations such as band limiting and matched filtering may be 
performed in this plane. The theory is developed and some applications 
presented in Refs. 3 and 4. 

By axis translation and successive transforming, other processes 
such as cros s -correlation , auto-correlation, convolution, spectral 
analysis, and pattern analysis are possible. Many of these are discussed 
in Refs. 3, 5, and 6. Specifically such systems have found application 
in such diverse areas as character recognition [3 and 7], image improve- 
ment [8], synthetic aperture radar [9], and antenna simulation [6]. 

C. SCOPE OF THE THESIS 

This thesis was conceived as a study of the merits of a low cost 
coherent optical processing system for student research. Although most 
systems in use perform two successive transforms, a single -transform 
system was considered more basic to the problem and such a system was 
constructed. Figure 2 is a photograph of the system. When completed. 
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Figure 2. The experimental system. The noted components are: (1) 1.5mw He-Ne 

laser and condenser lens in a light baffle; (2) pinhole; (3) collimating lens; 

(4) the input plane; (5) the transform lens; (6) shutter; (7) output plane with film 
pack. 




the system was used as a tool to investigate two otherwise unrelated ^ 
topics; These topics were a study of film grain characteristics and a 
proposed method of ship detection via optical detection of the ship wave. 

D. THESIS OUTLINE 

Section II considers the effects of film in a single transform optical 
system. A model of the film grain is proposed and a statistical distri- 
bution of grain size is assumed. From this, a semi-quantitative mathe- 
matical prediction of the grain transform is developed. The grain size 
distribution of Polaroid Type 55P/N film was determined with photomicro- 
graphs and a prediction of the grain' transform was made . The grain 
transform prediction was verified with an actual transform of a film 
specimen obtained with the constructed optical system. 

A method of detecting ship waves is considered in section III. A 
relationship between light reflected from a ship wave to the camera and 
the optical transform of the resulting photograph was developed sufficiently 
to verify that a deterministic relationship does exist. A sample ship wave 
photograph was then operated on by a single Fourier transform with matched 
filtering in the frequency plane. The matched filter was rotated to inves- 
tigate the-possibilities of further electrical filtering of the frequency plane 
, > 

output. 

Some suggestions for further research are given in section IV. 
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II. PROPERTIES OF FILM IN A SINGLE TRANSFORM SYSTEM 



A. INTRODUCTION 

Photographic film is the most commonly used medium for the input 
and output in an optical system. The primary advantage of film is its 
high information capacity which Jones [10] has compared with other media 
considered as information channels. The three film characteristics 
considered in this section may be thought of as disadvantages or limiting 
factors. The first factor is the time required to process the film, which 
in some cases may run several minutes both for the input and output. 

There are presently two means available to reduce this time; a Polaroid 
type film has a total processing time measured in seconds, while most 
conventional films can be processed in a specially prepared solution 
called a monobath which may require a few minutes. Based solely on 
time, the Polaroid type film seems the obvious choice. Although Polaroid 
negatives require clearing and emulsion hardening before use, the process 
could easily be mechanized and the time made very short. The Polaroid 
Type 55P/N film which provides both a positive and a negative was select- 
ed for evaluation in this project. 

The two other limiting factors are the film's grain structure and its 
non-linear recording characteristics. These two factors were investi- 
gated in particular for Type 5 5 (negative) film, but the method is quite 
general . 
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B. A FILM GRAIN MODEL 



Photographic films are generally constructed as shown in Figure 3. 

The process of exposing, developing, and fixing the film has the end 
result that the exposed grains are opaque and the unexposed grains have 
been washed from the film. This is described more fully by James and 
Higgins [11]. Even with no exposure, there are a certain minimum 
number of grains which will remain on the film as opaque spots. This, 
together with the density of the base, sets a minimum achievable density 
which is often referred to as fog level. The number of grains present 
increases as the exposure level increases and achieves a maximum den- 
sity at some level dependent upon the particular film. Since this grain 
will be Fourier transformed and imaged along with the signal, its effect 
is of some importance. 

The effects of grain have been treated in several papers such as 
Refs. 12, 13, and 14, where the grain is considered as noise having 
statistical properties. The problem is then one of a signal in the presence 
of noise and the techniques of communication theory are applicable. This 
treatment is ideally suited to double transform systems. However, this 
section is concerned with a single transform system, and the spatial 
frequency content of grain is of more interest. 

A review of literature disclosed no appropriate model of film grain 
for a single transform. It became necessary to develop a model and test 
its validity in a single transform system. This model is developed in the 
following paragraphs and though it did not result in a completely mathe- 
matical expression, the result is seen to be of value. 



Z2 




emulsion 



Figure 3 . 



Cross-section of a negative film. 
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Figure 4 is a photomicrograph of Type 55 (negative) film with no expo- 
sure. The density is thus D . which is specified 0.18, The photograph 

min 

suggested a grain structure model of overlapping circular stops. Since 
the diffraction pattern due to apertures distributed j identically to the 
stops will be essentially the same, the model is initially composed of 
apertures which is mathematically easier to handle. 

Francon [2] considered the Fourier transform of several apertures 
identical in shape. The model based on Figure 4 is composed of circular 
apertures of varying size but with the size-frequency distribution constant. 
A grain class was defined as a set of grains (or apertures) having one 
specified radius . The work of Francon may now be used to find the trans- 
form of each class. Thus: 

Na j Y" (^2^in + ysYin) 

fa(^2/y2)=f(^2/y2)a ^ © (2-1) 

n=l 

where Na is the total number of grains in class a, f(x 2 ,y 2 )a 1^ the Fourier 
transform of one aperture of class a, fa(>^ 2 /y 2 ) is the Fourier transform of 
the set called class a, k equals 27 t/X , and F is the transforming lens 
focal length. In words, the output due to the set of apertures in class a 
is the summation of Na outputs, each identical in magnitude but of varying 
phase. The total output of a given area of film becomes the summation over 
the area of all grain classes. 

^T 

C(X2,y2) = 2 fa(>^2/y2) (2-2) 

a=l 
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Figure 4, Photograph of Polaroid Type 55 
(negative) film with no exposure. (lOOOX). 
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where is the total number of classes, and f^(x 2 ,y 2 ) is the total Fourier 
transform. The intensity in the transform plane is: 



I = 






T j - (xsXin+ysyin)^ 

= I £ f(X 2 ,y 2 )a 2 : e 

a=l n=l 



(2-3) 



N, N, 



- r - «t - B J p r^2(^in ^im)"*‘y2 (yin"yim | 

r {f^(x2,y2)a2: r e ^ -^ | 

a=l n=l m=l 



For Na large, there are a large number of sine and cosine terms , the sum 
of which goes to zero except when m=h. Thus, this term sums to and. 



a, 

I = r [f® (Xj.yj. N.1 
^ a=l 



(2-4) 



But this is not the complete story. From equation (2-3), the intensity at 
the X 2 ,ys -origin iis ; 



I 

P 



a-p Na Na Siji 

= I r£=(x3,y2),r r ll= r |£^(x2,y2).N.= 

a=l h=l m=l a=l 



(2-5) 



The intensity at the origin is a function of the square of the number of 
grains rather than the number of grains as- in equation (2-4). Francon [2] 
indicates that, for equal apertures, experiment shows this bright central 
region to be surrounded by a dark : ring. Since there are a large number 
of apertures in a limited surface, there is necessarily a certain distri- 
bution of separation between each grain and its nearest neighbors which 
is not entirely random. The resulting order in the separations produces 
this dark ring. In the same reference, Francon shows that the remainder 
of the output is granular in appearance since the intensity deviation due 
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to the phase terms is of the same order as the mean. This effect is noted 
in many other papers and is the cause of the granular appearance which 

is characteristic of coherent illumination. 

The transform of a single circular aperture was then considered. This 
is covered in some detail by Goodman [3] and is: 



f(>^2/y2)a = e^LcircCo^a )] = O' a |_' 



r2Ji(z) 



(2-6) 



where Ca is the radius of apertures in class a, Ji(Z) is the first-order 
Bessel function, and fg is the incident light function (a constant). The 
other terms are defined as follows: 



circ(aa) 



1 for 

0 elsewhere 



Z = ^a.p (p“ = x/-'-y|) 

(This function has circular symmetry about the optical axis and is sketched 
in cross-section in Figure 5). The intensity distribution due to one 
aperture is the square of the function in Figure 5. Figure 6 is the trans- 
form of a 100 micron radius aperture and was obtained with the experiment- 
al system described later. From equation (2-5), the expression for 
intensity now becomes: 

( 2 - 7 , 

a=l 

This neglects the central region and the granularity effects . 

The original intent was to follow this line of analysis to find an 
expression for the total transform with the size-frequency distribution of 
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Figure 5. The first-order Bessel function. 
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Figure 6, Fourier transform of a 100 micron 
radius aperture, (focal length=50cm) , 
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a a known. However, this involves the distribution function of the out- 
put which is obtained by finding the inverse of the term 2 ]-^ (Z)/Z. No 
practical method of obtaining the inverse was found. The output can be 
obtained rather easily, though, for a finite number of classes by graphical 
summation. This was the technique used to complete the analysis. 

Up to this point, circular apertures have been considered, whereas 
the grain structure model involves circular stops . Defining fgs Qs the 
transform due to a circular stop model, fgj^ as the transform due to a 
circular aperture model, and fo as the transform of the incident light 
(source light), the following relationship applies: 

fss(^2/ys) = fo(^2/ys) - f2a(^2/y2) • (2-8) 

% 

The function fo(x 2 ,yg) will be the transform of a circular aperture defined 
by the system's limiting aperture. .This is generally the beam size itself, 
and since it is large compared to a grain size, the pattern fj, will be 
confined to the region near the origin. At a small distance from the 
origin, only fg^ and fgg are significant and: 

^2S (^2/y2) “ “ f2a(^2/y2) (2“9) 

or: 

^2S (^2^y2) ~ ^2a (^2/y2) (2 “10) 

That is to say that the intensity distribution of the transforms of comple- 
mentary inputs are identical except at the origin. This is Babinet's 
Complementary Aperture Theorem [2] . 

Considering all of the above factors, the film grain transform will 
exhibit a brilliant spot at the origin due to the function f^, (x 2 ,y 2 ) (in 
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effect - the DC term), and the factor of equation (2-5). Surrounding 
this is a dark annular ring as noted above. The remainder of the pattern 
will be granular in appearance with a mean distribution given by equation 
(2-7). 

Thus far this analysis has ignored the fact that the film thickness 
varies. Although some variations will occur solely because of manu- 
facturing tolerances, the most prominent effect is due to the processing 
action of exposed grain removal [3] . By defining , y^) to be the 

additional phase introduced at each grain center, equation (2-1) becomes: 



The analysis then follows as before with a similar result. The output 
pattern is essentially the same as given in equations (2-4), (2-5), and 



equal to zero, and some of the intensity at the origin must be dispersed 
over a slightly larger area. The dark ring surrounding the origin will not 
necessarily be at the same radius since the causal order depends on both 
the grain-center separations and the emulsion thickness variations. 

C. THE PREDICTED TRANSFORM RELATED TO FILM DENSITY 

A size-frequency distribution was obtained from Figure 4 and the 
results given in Table I and Figure 7. A size count such as this is some- 
what arbitrary and very dependent on the judgement of the observer. The 
guidelines used were: (1) for clumps of more than one grain, choose the 



Na 

fa(^2/y2)=f(>^a*y2)a ^ 

n=l 




(2-11) 



(2-7) except near the origin. Note that for Xg=y2=0 , (]) is not necessarily 
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Number of grains 




Figure 7. Size-frequency curve of 
Polaroid Type 55 (negative) grain. 



32 



TABLE I. GRAIN SIZE COUNT. 



Class 

ameter 


Diameter 

Range 


Number 


Weight (num 


0.4iU 


0.3-0.5/i 


259 


.234 


0.7 


o 

f 

o 

• 

00 


377 


.341 


1.0 


0. 9-1.1 


315 


.285 


1.3 


1 .2-1 .4 


84 


.076 


1.6 


1 .5-1.7 


47 


.042 


1 .9 


1.8-2 .0 


15 


.014 


2.2 


2. 1-2 .3 


4 


.004 


2.5 


2 .4-2 .6 


5 


.005 


2.8+ 


2.7- «> 


0 


.000 




Total 


1106 


1.001 



largest sub-units which could be considered as circular and count each 
as one grain, (2) whenever possible, choose sub-units such that overlap 
is minimized. The resulting curve of Figure 7 is similar to a typical 
curve given in Ref. 11 . 

The data of Table I was used to determine the intensity distribution 
in accordance with equation (2-7). Figure 8 is a plot of the transforms of 
each size-class and the total transform. There are no sharply defined 
nulls or peaks in the total transform which was to be expected since this 
is the transform of a continuously distributed input. The curve is a 
smoothly decreasing function of radius, which is difficult to study by 
experiment without some means of accurate light intensity measurement. 
Photographic film does not supply the desired accuracy and, at this time, 
no other means was available. At this point, it seems worth mentioning 
that the transforms of many signals of interest are not continuous spectra. 
As such their patterns will approach point spectra which are more easily 
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Figure 8. Fourier transform of Type 55 (negative) film grain 
by graphical summation." (focal length=10 cm). 
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analyzed. Film can provide only a rough quantitative analysis for the 
continuous distributions as in the present case. 

Type 55 (negative) film was chosen for the output medium. The 
D/logE curve (see Appendix A) is non-linear except in the central region 
of the slope. This is generally true of any film. In the linear region, 
the following relation holds: 

D = y log E + constant 

= y log IT + constant (2-12) 

where E is the photographic exposure, and T is the time duration of the 
exposure. This relation is valid from D=0.375 to 1.5. The range of 
output intensities far exceeded the linear range of the film and many 
exposures were needed at various values of T to cover the range. Thus, 
by considering the linear density regions of output films exposed at 
various values of T, the effective linear density range was extended. 
Since the intensity at a given radius must be a constant for all photo- 
graphs, two successive photographs were compared by letting: 

AD = y log 

t2 J- 2 

and set Ii=l 2 . Thus: 

AD = y log i (2-13) 

■ 1-2 

In practice, a AD was determined by comparing densities at equal radii 
in two photographs . 
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D. THE GRAIN TRANSFORM BY EXPERIMENT 



The optical system for this experiment is diagrammed in Figure 9. 

The position of and are reversed compared to the normal arrangement, 
which alters the transform kernal from exp(jk/F) yg) to expOk/d^) 

(xjXg+y^yg) and offers some control over the size of the output [3]. The 
input no-exposure transparency was placed against an optical flat with 
a matching fluid (mineral oil) between the two. Several exposures of the 
transform were taken to cover the intensity range with the input trans- 
parency base side against the fluid. A similar set of exposures was made 
with the emulsion side against the fluid. A single exposure was taken 
with the input transparency unmatched by the fluid. A visual comparison 
was made with the single unmatched output and the two matched cases at 
the same exposure level. The case with the fluid against the film base 
appeared virtually identical in density to the unmatched case, indicating 
that the base thickness variations have a negligible effect for this film. 

The emulsion-matched set of exposures correspond to equation (2-7) with 
no random phase variations. The base-matched set of exposures corre- 
spond to equation (2-11) because of the phase variations due to the 
unmatched emulsion side. 

The density of each of the exposures was recorded using a micro- 
densitometer to track a diametric line through the origin. Figures 10 and 
11 are typical density recordings from each of the sets. Figures 12 and 
13 are the corresponding Type 55 positives for these two recordings. 
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9. Optical system for grain structure Fourier transform. 
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Figure 10. Microdensitometer recording of Type 55 (negative) film grain transform with emulsion 
surface fluid matched. 
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Figure 11, Microdensitometer recording of Type 55 (negative) film grain transform with base 
surface fluid matched. 




Figure 12 , Transform of 
grain in Type 55 (negative) 
film with emulsion surface 
fluid matched. (Density 
recording is Figure 10). 



Figure 13. Transform of 
grain in Type 55 (negative) 
film with base surface 
fluid matched. (Density 
recording is Figure 11). 
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The exposure level is not the same for these two photographs and no 
comparison of the exact density levels should be made. They are pre- 
sented here only as typical of each set. In both figures, the peak is in 

saturation while the minimum level is below D , . The bright central 

min 

peak and dark annular ring are obvious in Figures 11 and 13. In Figure 
12, a faint ring is visible at a larger radius and appears as two spikes 
in Figure 10. Although this ring may also be due to the non-random 
order of apertures separated by short distances, it would be impossible 
to verify it without further knowledge of the film construction. A fine 
structure is evident in the recordings but the true character of granularity 
cannot by recorded by the microdensitometer since the variations occur 
in the order of one wavelength (0.6328 microns). The "tailing-off" in 
both figures is due to the film non-linearity in the toe of the D/logE 
curve . 

The series of exposures in each set were compared in density as 
noted above. The resulting density curves are plotted in Figure 14, The 
theoretical curve of Figure 8 was converted to density levels and also 
plotted in Figure 14. Note that the theoretical curve was derived from 
equation (2-7) which does not consider the effects near the origin. It 
must also be in error due to the effect of overlapping grains . When the 
grains overlap, the resulting shape is no longer circular and the dimension 
in one direction becomes much larger on the order of multiples of the 
smaller dimensions. Referring to Figure 8, note that the transforms of 
the larger grains affect the resulting curves at low values of radius . The 
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Figure 14. Normalized Fourier transform curves for Type 55 
(negative) film grain. Density levels based on a linear 
D/logE characteristic using Type 55 (negative) to record 
the output, (focal length=10cm) . 
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small grain transforms have a proportionately greater effect at large 
radii. The photomicrograph of Figure 4 shows very little overlapping for 
the smaller grains, keeping in mind the judgement criteria upon which 
the grain count was based. Thus, the experimental results at large 
radii should closely approximate theory. The curves of Figure 14 were 
normalized by equating the density levels at a radius of 12 cm. 

The experimental and theoretical curves show excellent agreement 
between 5 and 12 cm radius. There is no reason to believe that they 
would not continue to agree at larger radii but the curve could not be 
extended further because of lens curvature of field. The density departed 
from the theoretical curve for small radii in the manner predicted. That 
is, the intensity was higher near the origin due to the strong DC term 
and the effects of overlapping. It was noted also that the peak of the 
base-matched curve was less than for the emulsion-matched curve and 
that the intensity was somewhat more dispersed about the origin. 

The method of this section has been both quantitative and qualitative. 
Exact agreement between theory and experiment is not possible because 
the model is not an exact copy of the true grain structure . Had a better 
means of determining the intensity distribution been available, the 
normalizing procedures would not have been necessary and the results 
would, admittedly, carry more validity. However, the close agreement 
which does exist between theory and experiment is sufficiently convincing 
to prove the essential correctness of the model. 



43 



E. THE SIGNAL-PLUS-GRAIN TRANSPARENCY 



Thus far only the transform of the minimum grain level has been 

considered. This must now be related to a transparency composed of 

this minimum grain density with a signal superimposed on it. A typical 

transparency may consist of a continuous distribution of densities or a 

set of binary density levels, D and D . , but in either case the 

max min 

signal is on a background field of D . or D . Since there exists a 

min max 

limited density range, the signal itself must consist of some grain 
structure even in the binary case. 

Figures 15 and 16 are photomicrographs of Type 55 (negative) film 
with density levels greater than and they should be compared with 

Figure 4. It is obvious from Figure 14 that the circular grain model is far 
less representative of the actual grain structure when the grain structure 
is composed of irregular clumps. The nature of the exposure-development 
process is not fully understood, though several plausible theories have 
been proposed [11]. It is important here to note only that a clumping 
effect does occur and that it increases with density. Larger dimensions 
are now involved and, from the previous discussion in the effect of 
increased dimensions, an increase in intensity near the origin occurs at 
the cost of intensity at large radii. Assuming the input transparency 
area to be large compared to the size of any one grain, the size distri- 
bution is identical in all directions, and the output transform due to grain 
will still be symmetrical about the origin. The central DC term is reduced 
at the same time due to the overall increase in density. 
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Figure 15, Photomicrograph of Type 55 
(negative) film. Density 0.7, 



Figure 16. Photomicrograph of Type 55 
(negative) film . Density 1.7. 
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The shifting of energy to smaller radii due to. larger grain dimensions 
continues until the total grain area of the negative is approximately equal 
to the transparent area. Further increases in exposure produce negatives 
whose grain structure is approximately the complement of that for the 
negatives with lower densities. That is, the structure is one of trans- 
parent apertures rather than opaque stops. Indeed, in Figure 16, the 
structure appears to be a rough complement of Figure 4. The Complement- 
ary Aperture Theorem of Babinet predicts identical diffraction patterns for 
complementary screens except for the central DC term. Exact complements 
cannot be expected because the apertures of the high density negatives 

f 

do not have the circular structure of the stops of the low density nega- 
tives . For. these structures, an increase in density results in a shift of 
energy in the transform from small to large radii. The central DC term 
continues to decrease. 

A signal transparency will contain more than one density level and 
the output transform due to the grain alone will be dependent on the 
existing input density levels and relative areas which they occupy. The 
binary input consists of areas of complementary grain structures and the 
output will be the signal transform plus a grain transform which approaches 
that in Figure 8 or 12 . The worst case occurs when the input consists of 

.i 

'1 

intermediate densities where the grain size is about equal to the aperture 
(between grain) size. The transform intensities are more dispersed. 

Returning to Figure 14 and remembering that photographic film has 

a limited recording range (AD is about 1.5 for Type 55), it is seen that an 

adjustment of exposure level to place the tail of the curve below D . 

min 
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leaves only the grain transform near the origin. The effect of grain (with 
the emulsion matched) in the transform would be a central disc on the 
order of 0.2 cm radius for a focal length of 10 cm or, in general, 0.02X 
(focal length in cm) for Type 55 (negative) film. This means that the signal 
can have no significant frequency content less than about 2x10^ cycles/cm 
for this film . 

Figures 17, 18, and 19 are a series of signal transforms which were 

designed to test the above hypothesis and serve as a visual demonstration 

of transforms in the presence of grain. In practice, the best output was 

obtained at an exposure level which included some of the tail of the 

density curve of Figure 14. Figure 17 is the transform of grain 

alone. The signal input for Figures 18 and 19 were 2mm high characters 

with density D . against a field of D . The fine structure of con- 
min max 

centric rings seen in all three photographs is due to a small aperture 
mask required to select a particular character on the input transparency. 
Although the grain transform does degrade the image, the distinguishing 
features of the signal transforms are readily apparent. 

F. CONCLUSIONS 

Although the experimental grain transform differed significantly 
from the theoretical grain transform, the essential correctness of the 
grain model was verified. The overall effect of the grain appeared as a 
reduction in bandwidth by setting a lower limit of useful frequency response 
in cycles/cm. This effect appears to contradict intuition since grain is 
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Figure 17. Transform of 
film grain. Exposure is 
the same as for Figures 
18 and 19. (focal length 
is 50 cm) . 



Figure 18. Transform of 
the character '4' in the 
presence of film grain, 
(focal length is 50 cm). 
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Figure 19. Transform of the character '5' 
in the presence of film grain, (focal length 
is 50 cm) . 
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usually thought of in terms of limiting the fine detail, or the high spatial 
frequencies in an image. However, the following discussion shows that 
this is a special case and the results are somewhat misleading because 
a relatively grainy film was used. 

It must be noted that no part of the signal, including the fine detail, 
can have dimensions less than the diameter of the smallest grain. In 
most cases', the lowest spatial frequency of even the finest detail in the 
signal will fall in the central region where the grain transform is signi- 
ficant. The higher frequencies of the signal are visible here only because 
the signal transform intensity is significantly greater than the grain trans- 
form intensity in this frequency range. As seen earlier, the total inten- 
sity of the grain transform is a function of the total number of grains (or 
apertures in the case of input transparency. The grain 

transform intensity was reduced by lim.iting the area of the input trans- 
parency to include little more than the signal area as in Figures 17, 1,8, 
and 19. Retention of the low frequency signal spectrum, including that 
in the fine detail would require that the signal transform intensity be 
greater than the transform intensity over the range of the signal bandwidth. 
Extremely fine grain films are required to disperse the grain transform 

spectrum and "flatten" the curve of Figure 8. 

Consider a limiting case, shown in Figure 20, in which the grain 
transform is so dispersed that it appears essentially constant over the 
spatial frequency bandwidth of the signal. In a typical signal, the 
transform intensity would be a decreasing function of frequency. Above 
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Figure 20. Idealized signal and grain spectrums 
for a signal on fine grain film. 



Transform 




Figure 21. Idealized signal and grain spectrums 
for a signal on coarse grain film. 
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some high frequency limit, the signal spectrum intensity will be less than 
the grain spectrum intensity. Under these conditions, the grain structure 
is seen to have a greater effect on the high frequencies or the fine detail 
of the signal. Figure 21 illustrates a case similar to that of Figure 20, 
but one in which the grains are coarse. This case resembles the experi- 
mental results. In Figure 21, the grain spectrum intensity is greater 
than the signal spectrum intensity at both low and high frequencies and 
effectively creates a pass-band. The high frequency limit was not ob-wiv 
served in Figures 18 and 19 because the intensity range in the output 
exceeded the film density recording range. Further, since the grain 
transform intensity obscured most of the important lower signal frequen- 
cies in Figures 18 and 19, reconstruction of the signal by Fourier trans- 
form of the output transparency was found to be impossible. 

Thus, film grain in a coherent optical system has an effect similar to 

noise in electrical circuits. The film grain spectrum in the limiting case 

of Figure 20 appears similar to Gaussian noise in that its spectrum is flat 

over the range of interest. It appears that in many cases involving 

binary signals on fine grain film, this noise could be assumed Gaussian, 

at least as an approximation. In the discussion relating to Figures 15 

and 16, it was noted that the film grain spectrum is a function of the film 

density for densities other than D . and D . This would mean that 

min max 

there would necessarily be some correlation between grain noise and a 
signal comprising intermediate densities. The assumption of Gaussian 
noise could then be, at best, an approximation. 
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III. DETECTION OF SHIPS BY KELVIN WAVE 



A. INTRODUCTION 

The wave pattern produced by a pressure disturbance moving on 
the surface of the water was first considered by Lord Kelvin [15] in 1891. 
For a disturbance such as a ship moving on a straight course, the wave 
height is considerable inside an angle bounded by two horizontal rays at 
il9.5 degrees with respect to the ship's course line [16], Further, the 
angle is not a function of speed. Such a pattern, when considered 
against a random sea wave pattern, is unique and optical Fourier trans- 
forming techniques may be applied for its detection. 

The method used in this project was the single transform of a wake 
photograph with a matched filter in the first transform plane. The matched 
filter was rotated on the optical axis and the resulting transmitted light 
measured with a CdS cell to examine the possibility of further electrical 
filtering. A thorough investigation was not conducted since this was 
considered as a feasibility study and as a practical use-test of the optical 
system . 

B. THE WAVE PATTERN 

The free surface equation of ship waves is very difficult to derive 
when exact hydrodynamical theory is used. The solution of Lord Kelvin's 
as well as later work are approximations [16 and 17], Even then the 
results are quite involved and not completely germane to the problem of 
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detecting the waves. Ursell [16] has predicted and plotted contours for 
a ship’s wave near the critical angle, 0^1 9. 5 degrees. Figure 22 is one 
of these plots , 

Along the critical angle, the wave crests are inclined at an angle 
of 54.75 degrees to the track of the ship. Outside the critical angle, 
the wave height falls rapidly to zero. The horizontal angle of inclination 
between wave crests and the track decreases from 54.75 degrees to zero 
degrees for angles from less than Be to points on the track. However, 
the wave height also decreases inside although not as rapidly as out- 
side 0c. Thus there is some horizontal inclination angle between 54.75 
degrees and zero degrees at which the wave heights become insignificant. 
By way of comparison, the wave heights near the critical angle die out as 
an inverse function of the cube root of distance while, in the interior, 
they die out as an inverse function of the square root of distance [18]. 

The wave height and period are functions of the ship's speed and, from 
Figure 22, it can be seen that the period also decreases with angle. In 
summary, as the angle decreases from to zero, the wave height 
decreases, the wave period decreases, and the wave crests’ angle of 
inclination with respect to the track decreases (normal angle increases). 
These factors allow a qualitative prediction of the Fourier transform of 
such a pattern to be of the form in Figure 23. 6^ is the direction of wave 

travel (normal angle) for wave components near B^, and 6 is the direction 
of wave travel for wave components which exist near the track. B equals 

/A 

90-54,75 degrees or 35.25 degrees. 
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gx/V 



Figure 22, Contours of egual magnitude of the dimensionless 
surface elevation pV due to a concentrated force 

Pq moving with velocity V, for gr/V^=27Tr/X near 100. (Here 
r denotes the distance from the disturbance, X the wavelength 
of waves moving with phase velocity V.) For the sake of 
clarity only positive values of the elevation are shown. [After 
Ursell, Ref. 16]. 




Figure 23. Approximate Fourier transform 
of a ship' s wave . 
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C. PHOTOGRAPHIC PROPERTIES 



Stillwell [19] analyzed the relation between the light scattered from 
a sea surface and the optical density of a scene photograph. The method 
is outlined here with some changes which are necessary to the present 
problem . 

Let the ship's wave height pattern be described by H(x ,y)^ where x 
and y are taken to be on the sea surface. Assume the sea surface is 
illuminated by a completely clear or a completely overcast sky. The 
illumination at the camera becomes: 

M(t)) = KN(M)r(^.) (3-1) 

where K is a constant, L(Q the luminance of the sky in the direction 
reflected into the camera, and Ti^o) is the reflectivity at the incident 
angle . Figure 24 describes the geometry for a one-dimensional surface. 
For the two-dimensional surface, and 6 are sim.ple functions of 6 and 
O’ . The water reflectivity is a function of the incident angle and this is 
illustrated in Figure 25. Thus the camera illumination under these condi- 
tions is a function of the sea surface wave height structure, H(x,y). The 
camera illumination can now be described as: 

M(x,y) = KiH(x,y) (3-2) 

where is a constant. (Note: The illumination is also a function of 
the angle over the field covered by the photograph. The illumination is 
maximum directly away from the observer and is present in all sea photo- 
graphs [19]. A narrow field of view will be assumed so that this factor 
is eliminated.) Using equation (2-12), the recorded film density is; 

D(xi,yi) = y log I(xj^ , yj )T + Ks (constant) 
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(3-3) 




Figure 24, One-dimensional sea 
surface reflection geometry. 




Figure 25, Water reflectivity. In the visible 
region, water reflectivity is related to the 
angle from the normal by the above curve. In 
addition to the wavelength dependence of the 
reflectivity, the vertically (V) and horizontally 
(H) polarized components have different reflec- 
tivities, [After Stillwell, Ref, 19], 
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where T is the exposure time and the coordinates / are in the 

film plane. Figure 26 illustrates the relation between the surface co- 
ordinates, (x,y), and the film coordinates, (x^ , ) . The magnification, 

m, is a function of the film to object distance. Assuming r^+rg is 
approximately equal to T-y+Vz over the region recorded, then m may be 
assumed constant over the region and simple trigonometric relations may 
be applied. This approximation is generally true since the camera height 
is' usually much greater than the focal length. The coordinates , x and x^, 
are arbitrarily assumed to be in horizontal planes such that the relation, 
Xj^=mx, holds. The coordinates, y and y^ are then in a vertical plane 
containing both y and y^. The plane, P', is drawn normal to the line of 
sight and the coordinates related to this line become: 



mx 



yi = my 



and 0 < m < 1 



(3-4) 



The following angles may be derived. 



e = 90 - a 



c = a - sin 



, -1 yi/m 



r 



for r > > r 




r 



And finally, the Law of Sines, 






_y_ 

sin i> 



(3-5) 



Sin c 
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Figure 26, Sea surface-film plane geometry. 
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gives the required relation. 



yi 



my sin g 

1 + y/r 



(3-6) 



For a high altitude camera such as in a satellite, a is approximately 90 
degrees and equation (3-5) becomes. 



Yi = my 
and 



(3-7) 



= mx 



D. THE SHIP WAVE TRANSFORM 

The exposed photograph contains the sea surface information in the 
form of density variations as given by equation (3-3). This transparency 
is inserted into the input plane of a single Fourier transform system. The 
transmitted light amplitude is given by: 

, E),(^i / Yi) ■ 

fi(Xi,yi) = fo 10- ^ (3-8) 

The light intensity function in :the output or transform plane is: 

I (Xs/Ys) = ks |^[fi (Xi,yi)] (3-9) 

Combining equations (3-2, 3, 8 , and 9), the final expression is: 

I (x 2 ,Y 2 ) = K 4 l^[H‘(x,y)- V2 ] |2 (3_10) 

where: K 4 = K 3 (fo lO^^KiT)^ and x,y are 
given by equation (3-6) or (3-7) 

Although the above method is not complete and includes several assump- 
tions , the results are sufficient to verify the existence of a predictable 
relationship between a ship wave and the Fourier transform of its ,1 - 



photograph . 



60 



In a practical situation the ship wave must be considered in the 



presence of sea waves. The sea waves will not have the unique Fourier 
transform of the ship wave and, on the basis of this alone, should be 
easily separated in the transform plane. However, the sea waves do 
contribute to an overall noise level and it may be assumed that there is 
some minimum ratio of ship wave height to sea wave height below which 
the optical detection would be useless . Later mention will be made of 
this in section IV. 

E. A SHIP WAVE DETECTION EXPERIMENT 

The photograph of Figure 28 is a ship wave in a low sea state. The 
characteristic features are readily apparent along the critical angle. The 
transparency for Figure 2 8 was inserted in the input plane of the optical 
system illustrated in Figure 27, and an exposure was made in the transform 
plane. This is shown in Figure 29. The transform was similar to that 
predicted in section III-B. The differing density of the two pairs of wave 
lines may be traced to Figure 28 where the location of the sun produced 
differing modulation levels on the two sides of the wake. It was noted 
that this affected primarily the density rather than the shape of the out- 
put transform. The remainder of the pattern is composed of the grain 
transform and the central DC intensity considered in section II, and the 
transform of the sea waves which appear as noted in Ref. 19. 

A crude matched filter was produced by cutting out the portion of 
Figure 29 which represents the ship wave alone. This filter was then 
aligned in a mount which allowed rotation of the filter about a plane 
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Figure 28, Ship wave in 
a low sea state. 




Figure 29. Fourier transform 
of the ship wave of Figure 28. 
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normal to the optical axis. The mounted filter was placed at plane, P2 , 
and a CdS cell covered with diffusing glass placed behind it-. The geom- 
etry was designed so that the area of the transform pattern matched the CdS 
cell face. The cell output was then read as a function of angular filter 
position . 

This particular experimental arrangement was chosen to investigate 
the possibility of performing additional filtering by rotation of the matched 
filter at a constant rate. The CdS cell output would then be a time cor- 
relation of the signal and the matched filter. However, it must be noted 
that a CdS cell has a slow response time and could not be used at prac- 
tical rotation speeds. Referring to Figure 29, the output should consist 
of an average level plus components due to sea waves and the ship wave 
with a fundamental frequency of twice the rotation rate . The ship wave 
component would have a shorter time duration and produce different fre- 
quency harmonic terms. The resulting data is presented in Figure 30. 

This is the best run of a series. The data not included here was similar 
to that of Figure 30 but, in general, was degraded by inaccurate filter 
construction and filter centering. The construction problem could have 
been eased by using a photographic transparency filter, however, this 
would also have introduced undesirable film grain effects. The portion 
of the curve in Figure 30 representing the ship wave is quite prominent 
and suggests the possibility of visual detection. However, since the 
intended use of the output would involve frequency analysis, the curve 
was analyzed further for frequency content. 
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Figure 30. Ship wave detection output voltage curve. Each 
point is the average of V 0 and V 0^-2 8 O minimize the effects 
of misalignment. 







Figure 31. Ship wave detection. Fourier series coeffi- 
cients of output electrical signal. (Curve joining points 
is only for clarity.) 
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The FORTRAN program, FORIT, which computes Fourier series coef- 



ficients for a tabulated, periodic function was used for frequency analysis. 
Nineteen equally spaced points over 180 degrees of Figure 30 were used 

f' 

to obtain coefficients to the ninth harmonic. For comparison a similar 
curve, in which the ship wave component was replaced by a straight 



line joining the component end-points. The resulting exponential Fourier 



series coefficients are presented in Table II. The values normalized to 
the DC level 



TABLE II. EXPONENTIAL FOURIER SERIES 
COEFFICIENTS OF THE CdS CELL OUTPUT. 



HARMONIC 


COEFF. WITH 
WAKE 


COEFF. W/O 
WAKE 


NORMALIZED 
WITH WAKE 


NORMALIZED 
W/O WAKE 


0 


2 . 3341 V 


2 . 2552 v 


1.0000 


1.0000 


1 


0.0927 


0.2084 


0.0397 


0.0927 


2 


0.0871 


0.0883 


0.0373 


0.0397 


3 


0.1665 


0.0881 


0.0713 


0.0392 


4 


0.0995 


0.0690 


0.0426 


0.0307 


5 


0.0726 


0.0263 


0.0311 


0.0117 


6 


0.0679 


0.0392 


0.0291 


0.0174 


7 


0.0642 


0.0486 


0.0275 


0.0216 


8 


0.0490 


0.0038 


0.0210 


0.0017 


9 


0.0543 


0.0179 


0.0233 


0.0080 



are given in Table II and plotted in Figure 31 . A marked difference was 
noted particularly in the fundamental and the second and third harmonics . 
It was concluded that the output is susceptable to frequency filtering to 
indicate the presence of a ship wave in the input transparency. 



F. CONCLUSIONS 

Thus , the combined optical-electrical detection of a ship wave has 
been shown to be feasible, at least on a crude level. Since the transform 
of a ship's wake is always centered on the optical axis, the system is 
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not dependent on the wave position in the input field. The rotation of the 



matched filter eliminates the dependence on angular positioning. A phase 
reference derived from the filter rotation drive could conceivably be used 
to determine the ship's course with a 180 degree ambiguity by comparison 
of the phase of the ship wave component in the output signal to the phase 
reference. These factors suggest the use of such a system in an auto- 
matic alarm system to indicate the presence of a moving vessel in some 
monitored area of the ocean. 
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IV. SUGGESTIONS FOR FURTHER RESEARCH 



A. FILM CHARACTERISTICS 

The film grain model developed in section II allowed a semi- quanti- 
tative prediction of the grain transform. The ultimate limitation of this 
method lies in the grain model itself. However, the results would have 
been more useful with a closed mathematical expression for the total 
transform with the input size-frequency distribution known. As stated 
before, from equation (2-7), this involves the distribution function of the 
output which is obtained by finding the inverse of the term 2Ji(Z)/Z. The 
expression for the grain transform could be obtained in a closed form if 
this stumbling block were overcome. Since the input size-frequency 
distribution is often known for a type of film, the results could be applied 
almost directly to predict the effects of grain for any film type. 

B. SHIP WAVES 

The ship wave detection problem has been treated here on a fairly 
basic level with several approximations and assumptions. Since the 
concept does have obvious applications in surveillance of the seas, it 
seems worthwhile to define the problem more fully. Some of the problems 
involved are: (1) the variation of reflected light modulation with the sun 

and observer positions, (2) the relation between a ship wave and sea 
waves travelling at various relative surface angles , and (3) the signal 
to noise ratio. 
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C . SEA SURFACE ENERGY SPECTRA 



In Ref, 19, Stillwell describes the use of a single optical Fourier 
transform system to determine the energy spectrum of sea surface wave 
heights or wave slopes. Research into this area is of importance not 
only to oceanographers but is also vital to studies in radar scattering and 
acoustic propagation. 

A brief outline of Ref. 19 is given here. As noted in section III, the 
light arriving at a camera above a sea surface is a function of the wave 
slopes because of the incidence angle dependence of water reflectivity. 
The density of the exposed film is then determined by the sea surface 
characteristics and the problem geometry. The optical Fourier transform 
of such an exposed film can be expressed in terms of the film density 
and source light intensity. Since the transform output records the spatial 
frequency components of the input, it will also represent the sea surface 
spatial frequency components with a scale or magnitude factor involved. 

It is shown that the spatial frequency components of the sea surface are 

related to the slope energy spectrum, which is expressed as; 

D (k/ m , ) 

<l(Ki^) = C 10 y sec^ ip (4-1) 

where C is a constant, D(k/m, ip) is the output or transform photograph 
density function, k is the radian wave number, m is the magnification 
factor, ip is an angle between the observer-surface normal plane and the 
observer-zenith plane. 0 = 0 on a line directly away from the observer. 
Stillwell presents some results by plotting density contours of the trans- 
form output. A typical contour is shown in Figure 32. 
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Figure 32. Density contours showing the long wave- 
length spectrum and spectral peak for a light to moderate 
wind of long fetch. [After Stillwell, Ref. 19]. 
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Several attempts at analyzing sea surface photographs were made 
using the present optical system* The original photographs were 8X10 
inch transparencies, of which an area of 4^X4-^ inches was selected and 
reduced to 1/2 0th size. The area selected is shown in Figure 33 and the 
transform output in Figure 34. A microdensitometer was used to obtain 
the contours in Figure 35. The results were generally unsuccessful but 
not totally so. The contours do indicate some pattern which could only 
be attributed to the sea surface structure, and they are roughly similar 
to those in Figure 32. Most of the difficulties can be traced to the 
quality of the input. The original 8X10 inch transparency was grainy and 
the reduction of a photograph usually entails some loss of high frequency 
content. Both of these factors are apparent in Figures 34 and 35, and they 
are serious enough to make the reduced negative virtually useless. 

It is suggested that further work might be done in this area using 
sea surface photographs taken for the prime purpose of use as inputs to 
the optical Fourier transform system. 
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Figure 33, Sea surface photograph. 



Figure 34. Fourier 
transform of sea sur- 
face of Figure 33. 
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D=1.4, 1.2, 1.0, 0.8, 
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APPENDIX A. FILM CHARACTERISTICS 



The typical D/logE curve has the form of Figure 36, A standard 
approximation for the negative type film is shown in Figure 37, The result- 
ing equation which describes the curve of Figure 37 is: 

D = y [log E - log i] for E(D . ) < E < E(D ) 

min mQX 



D 



max 



for E ^ E (D 



max 



(A»l) 



D . for E ^ E (D . ) 

min min 

Often the spectral sensitivity curve, gamma, D , and D . are the 

max min 

only values given. Equation (A-l) must then be used to determine the 

density as a function of exposure. Spectral sensitivity curves relate 

radiant light intensity (or energy density) as a function of wavelength to 

reach a specified film density, usually 0.5. The exposure, E, is the 

incident luminous flux, U, times the duration of exposure. Luminous 

flux and radiant flux, R, are related by [20] : 

luminous efficiency ( 77 ) = ^ (A-2) 

K 

which is 150 lumens/watt for 6328A. The exposure is then: 



E= -UT =_riRT _ 

unit area unit area 

The constant, i, is evaluated using the radiant light intensity for D=0.5 

from the spectral response curve. The result is: 

D = y [logE - log E(D=0.5)] +0,5 (A-4) 

for E(D . ) < E < E(D ) 
min max 
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Figure 36. Typical film characteristic curves. 
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Figure 37, Approximation of the film 
characteristic curve. 




When exact light levels are not known or desired, the following form is 
useful . 

D = y log E + constant (A“5) 

The effect of reciprocity failure is not considered here. Reciprocity 
failure, occurs when exposures involve a high intensity and a very short 
exposure time or a low intensity and very long exposure time. The effect 
is a reduction of film sensitivity and the exposure time must be increased 
accordingly [11]. 
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APPENDIX B. LIGHT INTENSITY METER 



The light intensity meter described here was found useful in con- 
junction with the optical system. It can provide a rough measure of the 
intensity distribution at any point in the system and eliminate some oi the 
guess work to obtain a proper exposure. 

The circuit is diagrammed in Figure 38 and is basically a high gain 
operational amplifier with a cadmium sulfide photocell in the feedback 
loop [23], This configuration was used to reduce the effects of temper- 
ature and power supply variations. Transistors and Qg act as a 
differential amplifier with a 1.5 volt input. Transistor Q 3 is a DC level 
shifter and amplifier while transistor Q 4 drives the meter. 

The photocell was mounted in a telescopic tube behind a 0.5mm 
diameter pinhole. The telescopic tube was affixed to a bracket with 
vertical and horizontal tracks so that the photocell position could be 
adjusted over the face of the bracket. The bracket slides into the film 
pack holder with the pinhole in the film plane. The completed unit is 
shown in Figure 39. 

The unit was calibrated from zero to 2990 meter-candles using 
diffuse illumination and a Weston Illumination Meter, There was no 
standard available to calibrate the meter above 2990 meter-candles. The 
calibration above 2990 meter-candles was accomplished by comparing the 
readings at two distances from the source. Letting lo (d) define a calibrated 
reading at distance d, the relation is: 
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Rl-500 ohm pot 

R2-2500 ohm pot 

R3-25,000 ohm pot 

R4-100,000 ohm pot 

R5-1 0,000 ohm, 

R673OOO ohm, tW 

R7-1000 ohm, 

R8, Rll- 470 ohm, ^ W 

R9-10,000 ohm, t W 

RlO-2700 ohm, \ W 

R -Clairex CL505-L 
pc 



Cl-lOuf , 15V elec . 

51- DPST, normally open 

52- SP, 5 pos . selector switch 
Bl-1 .5V 

B2-9V 

Ql, Q2-2N3391 

Q3-2N4126 

Q4-2N3391 



Figure 38, Light Intensity Meter circuit. 
[After Barnum, Ref. 21], 
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Figure 39. Light Intensity Meter 
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I. (di) 



constant, K 



(B-1) 



The distances , and dg , were used for all readings and the readings 
outside the calibrated range become: 



I Tu . ^ (di) = KIcCds) 
uncalibrated ^ 



(B-2) 



The resulting accuracy cannot be as close as ±5%, The calibration curves 
are given in Figure 40. A fourth meter scale was not calibrated because 
a light source of such high intensity was not available. This scale is 
beyond the range of intensities in the optical transform system and was 
not needed. 

In use, the probe is positioned to monitor the light intensity in a 
region of the output plane where a linear density response is desired. An 
approximately correct first exposure can be obtained by exposing for a 
density of 1.0 for most films [11], Equation (A-4) , solved for E to reach 



a density of 1.0, becomes: 



E 



E(D=0.5) 10^*^^^ 



(B-3) 
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Meter Reading 



Figure 40 . Light Intensity Meter calibration curves . Scale 
readings are: XI for Scale A(±5%), XIO for Scale B(±5% for 
meter reading less than 7 and 10% above 7), and XI 00 for 
Scale C (±2 0%). Functional tests: BATTERY-1 3 , 3 meter 
reading; CHECK-9.2 meter reading. 
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